Weak pseudogap behavior in the underdoped cuprate superconductors 
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We report on an exact solution of the nearly antiferromagnetic Fermi liquid spin fermion model 
in the limit nT 3> w s f, which demonstrates that the broad high energy features found in ARPES 
measurements of the spectral density of the underdoped cuprate superconductors are determined by 
strong antiferromagnetic (AF) correlations and precursor effects of an SDW state. We show that the 
onset temperature, T cr , of weak pseudo-gap (pseudoscaling) behavior is determined by the strength, 
£, of the AF correlations, and obtain the generic changes in low frequency magnetic behavior seen 
in NMR experiments with £(T cr ) ~ 2, confirming the Barzykin and Pines crossover criterion. 
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Magnetically underdoped cuprates may be distin- 
guished from their overdoped counterparts by the pres- 
ence of a maximum at a temperature T cr > T c in the tem- 
perature dependent uniform susceptibility, Xo{T). They 
are characterized by the occurrence of a quasiparticle 
pseudogap observed by NMR and INS experiments, op- 
tical, transport and specific heat measurements, and in 
angular resolved photoemission spectroscopy (ARPES). 
Barzykin and Pines [0 proposed that at T cr , sizable 
AF correlations between the planar quasiparticles bring 
about a change in the spin dynamics, and that at T 
near T cr , the quantity ^T^/ 63 ^, where ^ is the 
63 Cu spin-lattice relaxation time and 63 T2G is the spin- 
echo decay time, changes from being nearly independent 
of temperature (above T cr ), to a quantity which varies 
as (a + 6T)" 1 . This behavior has recently been con- 
firmed in NMR measurements by Curro et al. 0. Since 
63 TiT/ 63 TJq oc £" -z , where z is a dynamical exponent 
0, the near temperature independence of 63 TiT/ 63 T 2 G 
found between T cr and a lower crossover temperature, T„, 
suggests that one is in a pseudoscaling regime (z m 1), 
while the temperature independence of of &3 TiT /^T^q 
above T cr suggests mean field behavior (z sa 2). More- 
over, above T* ARPES experiments show that the spec- 
tral density of quasiparticles located near (71", 0) has devel- 
oped a high energy feature || . We refer to this behavior 
as weak pseudogap behavior, to distinguish it from the 
strong pseudogap behavior found below T*, where exper- 
iment shows a leading edge gap develops in the quasi- 
particle spectrum || . Strong pseudogap behavior is also 
seen in specific heat, d.c. transport, and optical experi- 
ments, while 63 T 2 g measurements show that the AF spin 
correlations become frozen (i.e., £ ~ const.) and e3 TiT 
displays gap-like behavior. 

The nearly antiferromagnetic Fermi liquid (NAFL) 
model ||,[5| of the cuprates offers a possible explanation 
for the observed weak and strong pseudogap behavior. In 
this model, changes in quasiparticle behavior both reflect 
and bring about the measured changes in spin dynam- 



ics. The highly anisotropic effective planar quasiparti- 
cle interaction mirrors the dynamical spin susceptibility, 
peaked near Q = (ir.ir), introduced by Millis, Monien, 
and Pines Q to explain NMR experiments: 
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V eS (q^)=g Xq M = 1 + g2(q _ Q)2 _.^ , (1) 

where xq = ct£, 2 , with a constant, and g is the coupling 
constant. 

Since the dynamical spin susceptibility Xq( CJ ) peaks 
at wave vectors close to (it, ir), two different kinds of 
quasiparticles emerge: hot quasiparticles, located close 
to those momentum points on the Fermi surface which 
can be connected by Q, feel the full effects of the interac- 
tion of EqJI]; cold quasiparticles, located not far from the 
diagonals, | fc^ | = \k y \, feel a "normal" interaction. Their 
distinct lifetimes can be inferred from transport experi- 
ments, where a detailed analysis shows that the behavior 
of the hot quasiparticles is highly anomalous, while cold 
quasiparticles may be characterized as a strongly coupled 
Landau Fermi Liquid [Q. 

In the present communication we focus our attention 
on temperatures T > T*. Our reason for doing so is 
that for T > T*, fits to NMR experiments show that 
c<j s f < ttT, the frequency equivalent of temperature. Be- 
cause of its comparatively low characteristic energy the 
spin system is thermally excited and behaves quasistati- 
cally; the hot quasiparticles see a spin system which acts 
like a static deformation potential, a behavior which is 
no longer found below where the lowest scale is the 
temperature itself and the quantum nature of the spin 
degrees of freedom is essential. We find that above T*, 
in the limit ttT 3> lu s { it is possible to obtain an exact 
solution of the spin fermion model with the effective in- 
teraction, V^ AFh (q,uj), of Eq. § 

Our main results are the appearance in the hot 
quasiparticle spectrum of the high energy features seen 
in ARPES, a maximum in Xo{T) and a crossover in 
63 TiT/ 63 T 2 2 G for £ > £ « v F /A a . These are produced 
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by the emergence of an SDW-like state, as proposed by 
Chubukov et al. || Here, vf is the Fermi velocity and 
A Q = \J (S 2 ) ~ | , a characteristic energy scale of the 
SDW-like pseudogap. Using typical values for the hop- 
ping elements (see below), and g ~ 0.6 eV (determined 
from the analysis of transport experiments in slightly un- 
derdoped materials 0), we find £ ~ 2. For £ > £ Q , the 
hot quasiparticle spectral density takes a two-peak form 
which reflects the emerging spin density wave state, while 
the MMP interaction generates naturally the distinct be- 
havior of hot and cold quasiparticle states seen in ARPES 
experiments. 

Before discussing these results, we summarize our cal- 
culations briefly. Using the effective interaction, Eq. [j], 
the direct spin-spin coupling is eliminated via a Hubbard- 
Stratonovich transformation, introducing a collective 
spin field S q (r). After integrating out the fermionic de- 
grees of freedom, the single particle Green's function can 
be written as 

G k Ar - r') = (G k , kCTCT (r, r'|S)) , (2) 

where G\^^ aa (r 1 r'|S) is the matrix element of 

[G£(t - r')5 k ,v ~ JLs k ^(r)5(T ~ r') ■ a}-\ (3) 

which describes the propagation of an electron for a given 
configuration S of the spin field, a is the Pauli matrix 
vector and G~£ — ~(d T + £k) is the inverse of the unper- 
turbed single particle Green's function with bare disper- 
sion 

£k = — 2<(cos k x + cos k y ) — At' cos k x cos k y — /x . (4) 

In the following we use t = 0.25 eV and t' = — OAt for 
the nearest and next nearest neighbor hopping integrals, 
respectively, and we adjust the chemical potential fi to 
maintain the constant hole concentration at = 15%. 
The average (• • -) = J VS ■ ■ ■ cxp {— S Q } is performed 
with respect to the action of the collective spin degrees 
of freedom: 

S o(S) = ^Xq'KlSqK) -S_ q (-iw„) (5) 

q,n 

where uj„ = 2mrT and Zb is defined via (1) = 1. In us- 
ing Eq. H we have assumed that (i) Xq(uj) is the fully 
renormalized spin-susceptibility taken from the exper- 
iment and (ii) any nonlinear (higher order in S than 
quadratic) terms of the spin field can be neglected. The 
model which results from assumption (ii) is usually re- 
ferred to as the spin fermion model. 

After inversion of Eq. |3] in spin space, the average of 
Eq. can be evaluated diagrammatically using Wick's 
theorem for the spin field. In the above mentioned static 
limit, irT w s f, it suffices to consider only the zeroth 



bosonic Matsubara frequency in x q (iu; n ). The remain- 
ing momentum summations are evaluated by expanding 
£k+q ~ £k+Q + V k+Q ■ (q — Q) for momentum trans- 
fers close to Q, using vg+Q = de^+Q/ 'dk a . In this limit 
all diagrams can be summed up by generalizing a solu- 
tion for a one dimensional charge density wave system 
obtained by Sadovskii || to the case of two dimensions, 
and more importantly, to isotropic spin fluctuations. We 
find the following recursion relation for the Green's func- 
tion Gk(w) = G k =0 (w), whose imaginary part determines 
the spectral density A(k,u>), seen in ARPES p0[ : 

[GLH]- 1 = w - e k+lQ ^A^H- (6) 

Here, i>k,z = |vk+cj| and m = (I + 2)/3 if I is odd, 
while Uk.z = |vk+Q|(|cos0 k | + |sin</> k |) and K[ = 1/3 
if I is even. 0k is the angle between Vk+Q and Vk- 
The recursion relation, Eq. @, enables us to calculate 
A(k, ui) to arbitrary order in the coupling constant g. 
In the limit £ — > oo the Green's function reduces to 
GkM = JdA p(A) G k DW (A), where Gjp w (A) is the 
single particle Green's function of the mean field SDW 
state and p(A) ~ A 2 exp(— |A 2 /A§) is the distribu- 
tion function of a fluctuating SDW gap, centered around 

Ao, i.e., the amplitude fluctuations of the spins S are 

confined to a region around |(S 2 ), although in our cal- 
culations directional fluctuations are fully isotropic and 
spin rotation invariance is maintained. Below we show 
that the SDW like solution is obtained even at finite val- 
ues of £. 

The quantities we calculate are the single particle spec- 
tral density, A(k, u>) and the low frequency behavior of 
the irreducible spin susceptibility x q (w,T). Our princi- 
pal results are depicted in Figures 1-3. The ^-dependence 
of the Fermi surface, shown in the inset to Fig. la, is 
similar to the results obtained by Chubukov et al. [|ll| at 
T = 0; however for the experimentally relevant range of 
£, hole pockets around k = (ir/2,ir/2) do not form. We 
have used the criterion £k + Re£k(<^ = 0) = to deter- 
mine the Fermi surface. Although strictly valid only for 
T = 0, for finite temperatures it indicates when a quasi- 
particle crosses the chemical potential. In the limit of 
very large correlation length, we find, in addition to the 
two broadened poles of a SDW like state, a third solution 
of lo = ReEk(w) + £k> which although not visible due to 
the large scattering rates, ensures that even for £ — > oo a 
large Fermi surface and only a pseudogap in the density 
of states occur. 

A comparatively sharp transition between the behav- 
ior of hot quasiparticles (located at points a and b on the 
Fermi surface of Fig. la) and cold quasiparticles (at d 
and e) is found. For hot quasiparticles the single par- 
ticle spectral density evolves with temperature as the 
AF correlation length increases from £ w 1 to 5; a two 
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peak structure, which corresponds to a transfer of spec- 
tral weight from low frequencies to frequencies above 
200 — 300 meV, develops at £ Q w 2, and is quite pro- 
nounced for £ > 3 (Fig. lb). As may be seen in Fig. 
la, the shift in spectral density found for hot quasiparti- 
cles does not occur for cold quasiparticles, whose spectral 
density continues to be peaked at the Fermi energy. 
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FIG. 1. (a) Spectral density multiplied with Fermi function 
on the Fermi surface for f = 3. The distinct behavior of hot 
and cold quasiparticles is visible. The inset shows the cor- 
responding Fermi surface and its ^-dependence. The dashed, 
solid and dashed-dotted lines correspond to £ = 1, 3 and 
10, respectively, (b) ^-dependence of the spectral density of 
point b in (a). The transition from a conventional line shape 
for small £ = 1 (dotted line) to spectral densities with broad 
lines and pronounced high energy features (£ = 2, 3, and 5 
for the dashed, dashed dotted, and solid line, respectively), 
reflecting a precursor SDW state, is visible. Note, the satu- 
ration of the precursor effects for £ « 3. 



vertex is considerably enhanced for the hot quasiparti- 
cles; it is almost unaffected for the cold quasiparticles, 
reflecting again their qualitatively different behavior. 

We combine the results for Gk(w) and the electron- 
spin fluctuation vertex function and so determine the ir- 
reducible spin susceptibility x q (w). We find (Fig. 2) that 
both Xo(T) and Xq(T) exhibit maxima at temperatures 
close to T cr where £ « 2. In these calculations we as- 
sumed that £ -1 (T) = \ + \ T ^ r ~_ T r. between T cr = 470 K 
and T* = 220 K and £" 2 (T) = \ + \ ^7?" above T cr 



\ tc7_t, between T 
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consistent with the NMR results of Curro et al. H for 
YBa2Cu408- The behavior of Xo{T) and Xq{T) above 
the maximum reflects the increasing importance of life- 
time (strong coupling) effects which act to reduce both ir- 
reducible susceptibilities. Because of comparatively short 
correlation lengths (£ < 2), it is likely that Eliashberg cal- 
culations j^] will provide a better quantitative account in 
this mean field regime. The fall-off in x (T) and Xq(T) 
below T cr arises primarily from the transfer of quasipar- 
ticle spectral weight to higher energies. 

The determination of the full spin susceptibility 
Xq(w) = Xq( w )/(1 — ^qXq( w )) requires calculating the 
restoring force J q , and is beyond the scope of the present 
work, since J q is determined by the renormalization of 
the spin exchange fermion-fermion interaction through 
high energy excitations in all other channels. However, 
for Xo(T), with the assumption that J q =oXo(2c r ) = 0.5, a 
good quantitative fit to the experimental results of Curro 
et al. between T cr and T* is found. 
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FIG. 2. Irreducible part Xi(T) of the static spin suscepti- 
bility at q = (0, 0) and q = (ir, n) as function of temperature. 



We sum all diagrams of the perturbation series for the 
electron-spin fluctuation vertex function in similar fash- 
ion as the Green's function Gk(w) in Eq. (^). The lack of 
symmetry breaking is essential for a proper evaluation of 
the vertex which, as long as the spin rotation invariance is 
intact, is reduced at most by a factor w | for the high en- 
ergy features juLO. For lower excitation energies, this 



We find that both vertex corrections and quasiparticle 
spectral weight transfer play a significant role in deter- 
mining the low frequency spin dynamics. As may be 
seen in the inset to Fig. 3, when both effects are taken 
into account, our calculated values of the spin damping 
7q = x'q(ui,T)/uj\^_ display the crossover at T cr an- 
ticipated by Monthoux and Pines |13[. A second calcu- 
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lable quantity which can be compared with experiment 
is Xq(^) 2 /7q = w s /xq, being proportional to the prod- 
uct, 63 TiT/( 63 T 2 g) 2 . As may seen in Fig. 3, qualita- 
tive agreement with the results of Curro et al. [p] for 
YBa 2 Cu40g is found. 

Physically, the most interesting aspect of our results is 
the appearance of SDW precursor phenomena, brought 
about by the strong interaction between the planar quasi- 
particles, for moderate AF correlation lengths, £ > £ w 
2, in contrast to earlier calculations, in which SDW pre- 
cursor behavior was only found in the limit of very large 
correlation length |l4|]. Our exact solution of the static 
problem enables us to access this region of strong cou- 
pling. For £ > £ , the hot electron mean free path, ~ 
begins to be small compared to £, so that the quasi par- 
ticle can no longer distinguish the actual situation from 
that of a SDW state; hence we find pseudo-SDW behav- 
ior, i.e. SDW behavior without symmetry breaking. The 
related shift of spectral weight for states close to (tt, 0) 
affects mostly the low frequency part of the irreducible 
spin susceptibility and leads to the calculated crossover 
behavior. Because symmetry is not broken, there will 
continue to be coherent states at the Fermi energy; al- 
though their spectral weight is small one still has a large 
Fermi surface. Such coherent quasiparticles, invisible in 
the present temperature range, are, we believe, responsi- 
ble for the sharp peak for k ~ (n, 0), observed in ARPES 
below the superconducting transition temperature H. 
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FIG. 3. Xq^/XqI _ as function of temperature, com- 
pared with experimental results of Ref. [2] for TiT/T^q. The 
inset shows the crossover in the calculated T-dependence of 

The present theory cannot, of course, explain the lead- 
ing edge pseudogap found below T», since one is then 
no longer in the quasistatic limit for which our calcu- 
lations apply Strong pseudogap behavior corresponds 
to a further redistribution of quasiparticle states lying 



within sa 30meV of the Fermi energy No appreciable 
change is seen in the high energy features found in the 
present calculations. It is likely that strong scattering 
in the particle-particle channel plays an increasingly im- 
portant role below T*, since we find above T* important 
prerequisites for its appearance; an enhanced spin fluc- 
tuation vertex and a pronounced flattening of the low 
energy part of the hot quasiparticle states In addi- 
tion, below T* the quantum behavior of spin excitations 
becomes increasingly important. This reduces the phase 
space for quasiparticle scattering, and leads to the siz- 
able suppression of the hot quasiparticle scattering rate 
found below T* @. 
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